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Abstract. It is shown how to extend the formal variational calculus in order to incorporate 
integrals of divergences into it. Such a generalization permits to study nontrivial boundary 
problems in field theory on the base of canonical formalism. 

1. Introduction. The Hamiltonian formulation of classical mechanics [Arn] is based 
on geometrical constructions which use such notions as differential forms, vector fields 
and multivectors, and such operations as differential, interior product. Lie derivative, 
Schouten-Nijenhuis bracket. Most of these constructions were extended to field theory 
in the process of studying nonlinear integrable models during the last 20 years [01v86]. 
This approach has been called the formal variational calculus [GD] because it ignores 
any terms arising as a result of integration by parts. This is fully justified in case of 
periodic boundary conditions or fast decay of fields at spatial infinity, but unfortunately, 
this method is not applicable in its initial form to many other problems interesting from 
physical point of view. For example, massless fields are slowly decaying at infinity and, 
as a result, some important characteristics of these fields are expressed just through 
surface integrals (or volume integrals of spatial divergences). They are necessary to form 
canonical generators of the global gauge transformations or asymptotic symmetries of 
the Riemannian metric. The great efforts were started at the end of fifties to understand 
the role of surface terms in the Hamiltonian of General Relativity [ADM]. Only after 15 
years of study the satisfactory explanation had been given [RT] . But even then not all 
questions were answered. For example, one might worry how to retain the surface terms 
which are necessary to realize the Poincare algebra in asymptotically flat space [RT], 
[Sol85] 

{H{0,H{T^)}^H{[tv]) 
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when doing local calculations of the constraints algebra 



{'H{x),H{y)} = g''\x)nk{x)5,a{x,y)-g-\y)Hb{y)^a{y,x), 

{H{x),Ha{v)} = -n{v)5,a{v.x), 
{Ha{x) ,Ub{v)} = Hb{x)d^a{x,y) - Haiy)S,biy,x). 



We will show in the following that all the main structures of the formal variational 
calculus can be extended to include nontrivial contributions of divergences through in- 
troduction of a new grading and new pairing compatible with it. So, it occurs possible 
to preserve the nice geometrical language in the more general case than before. After all 
the field theory Poisson bracket is given by a new formula which differs from the stan- 
dard one by surface terms. Simultaneously wc get the answer to the mentioned problem 
of disappearance of the surface contributions in local calculations with (5-function. The 
natural way to take the boundary terms into account is to introduce the characteristic 
function dn{x) of the integration domain il. Then relations like 



give the solution. In its turn this is connected with the observation [Sol92] that transfor- 
mations of the type (for example, transformation to Ashtekar's variables) 



in field theory are canonical only up to boundary contributions, because the standard 
Euler-Lagrangc variational derivatives in general do not commute [And76] , [And 78] . 

We expect that boundary conditions should be treated in this formalism as a kind of 
constraints put on the initial data, i.e., they should be added to the Hamiltonian with 
some Lagrange multipliers and then checked for compatibility with the dynamics. The 
requirement of compatibility may lead to secondary boundary conditions or to fixing the 
Lagrange multipliers. But now this subject is not enough studied and our consideration 
is preliminary and limited to one example: the nonlinear Schrodinger equation. 

2. New Poisson bracket formula. Below we use the local coordinate language and 

instead of the manifold with a boundary consider a domain D, in R" having a smooth 
boundary dQ. We do not expect that global formulation could meet with serious difficul- 
ties. 

Definition 1. An integral over a finite domain f2 of a function of field variables 
(t)"^{x), A = 1, ...,p and their partial derivatives Djcj)-^ up to some finite order 



is called a local functional. 

Remark 1. In contrast to the standard definition we do not treat these local 
functionals as equivalent if they differ by a divergence term. 

All the functions / and 0a as well as their variations throughout the paper are 
supposed infinitely smooth, i.e. C°°(R"). We use multi-index notations J = (ji, ...,jn) 




Pa{x) -^pa{x) + 



SFjq] 
6q^{x) 








Dj 



di^x^...di-x' 



.n 



Do = 1. 
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The derivative operator D will denote below the full partial derivative taking into account 
also coordinate dependence of fields (pA {x) ■ As the number of sums in some formulae of 
this paper is large enough we will write only a sign of summing without displaying the 
indices of summation. According to this rule, sum over all repeated indices should be 
understood. In those cases, where it is not so, we display the summation indices. Also, 
we do not show the limits of summation, because they are natural, i.e. outside them the 
summand is simply zero. Usually we omit d^x in the integrals and show the arguments 
only when they can be mixed. 

We denote as A the space of local functionals . It is important that this space in- 
cludes functionals with integrands depending on derivatives of arbitrary order [And92]. 
Otherwise the Poisson brackets could go out of A. The following is the general definition 
of field theory Poisson bracket. 

Definition 2 A bilinear operation {•, •} such that for any F,G,H e A 

1) {F,G}eA; 

2) {F,G} = -{G,F}; 

3) {{F, G}, H} + {{H, F}, G} + {{G, H}, F} = 0; 
is called the field theory Poisson bracket. 

The key idea of the new formula is in exploitation of the full variations which are free 
on the boundary. The variation of a local functional 

contains the differential operator which is called Prechet derivative 




We propose to define field theory Poisson bracket by the formula 

{F,G}=Y, I MDfjABDgJ, 
where trace of two differential operators is used 

Tr{AB)=Y,DiAKDKBi, 

I,K 

where 

A = J2AkDk, B = Y,BiDi. 

K I 

The important property of the trace is 

Tx{AD]3) = DTr(i^) = Tr(£)ii3). 

The individual structure of a Poisson bracket is given by matrix Iab ■ More general treat- 
ment of it will be given in the next Section and here it is simply constant antisymmetric 
matrix. 

Symmetrized covariant derivatives can also be used in the expression for the first 
variation of local functional . We can replace partial derivatives by covariant ones in the 
trace calculation if the curvature is zero or if one operator is simply a multiplication by 
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a function. Then covariance of the new formula under changes of independent variables 
is evident. In general case a special consideration is necessary. 

The new bracket differs from the standard one in the exact fulfilment of the Jacobi 
identity under arbitrary boundary conditions [Sol93]. In the same time its calculation 
is not more complicated than usual since we need not integrate by parts to get Euler- 
Lagrange derivatives. 

We can also use another representation of the first variation [01v86] 



where the higher Eulerian operators [KMGZ], [Aid] 

i^i(/) = E(-i)'"'^'"(j)^- 



K ^ ' ^ '-''Pa 

are used. The zero order operator is just the standard Euler-Lagrange variational deriva- 
tive. Binomial coefficients for multi-indices are 

where ordinary binomial coefficients are 



A = (ji/{k\{j-ky.) if o<A-<j; 

kj lo otherwise. 



Let us mention that if J is not contained in K, then all quantities having multi-index 
{K — J) are zero. The sums over J and K above are really finite because local functional 
can depend only on a finite number of derivatives according to Definition 1. 
A remarkable property of these operators is 

It leads to the result that if 



then 



{//,/.} ^ 
{J Df,J g} = {J f,J Dg} = J Dh. 



3. Extension of the formal variational calculus. In dealing with terms arising in 
the integration by parts it is suitable to represent integrals over finite domain as integrals 
over infinite space with the help of the characteristic function 

if X e n: 
otherwise. 

We can understand it also as Heaviside function 6q{x) = 9{Pq), where equation Pn{x) = 
defines the boundary and 

> if a; e O; 
< otherwise. 

Then 



I'-i 



enf, 
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and, for example, the full variation of a local functional can be expressed in the form 

SF ■ 

where the distribution 



1-0(PA-, 



^ = E\{enf) = Y,{-l)^'^DjenEi{f) 
d(j)A ^ 

could be called the full variational derivative. Such a representation corresponds to the 
situation opposite to the standard one: here distributions are of finite support whereas 
test functions 8<f)A are arbitrary. 

A grading in linear space L is a decomposition of it into direct sum of subspaces, with 
a special value of some function p (grading function) assigned to all the elements of any 
subspace [Dorf]. Elements of each subspace are called homogeneous. 

In our case the factor DjOn is responsible for the grading and the function p takes 
its values in the set of all positive multi- indices J = (ji, • • ■ , jn) 

j=o 

We always can return to the standard formal variational calculus by putting 9q.{x) = 1. 

A bilinear operation x,y .to y. defined on L, is said to be compatible with the 
grading if the product of any homogeneous elements is also homogeneous, and if 

p{xoy) =p{x) +p{y). 

^.1. Local Junctionals and evolutionary vector fields. Here we will call the expression 
given in Definition 1 the canonical form of a local functional. We formally extend that 
definition by allowing local functional to be written as follows 



oo 

E 

7=0 



Djen{,x)f('\ct>A{^),DKM^))d''x = Y, j O^'^ f^'\ 



where in accordance with the previous definition only a finite number of terms is allowed. 
Here and below we simplify the notation for derivatives of and remove il. Of course, 
any such functional can be transformed to the form used above through integration by 
parts with 

/ = ^(-l)l^lD./<^). 

So, the formal integration by parts over infinite space R" evidently changes the grading. It 
will be clear below that the general situation is from one side compatibility of all bilinear 
operations with the grading and from the other side compatibility of them with formal 
integration by parts. So, basic objects (local functional etc.) are defined as equivalence 
classes modulo formal divergences (i.e., divergences of expressions containing ^-factors) 
and the unique decomposition into homogeneous subspaces with fixed grading function 
can be made only for representatives of these classes. But we will see that the pairing 
will be defined in such a way to avoid any ambiguity. 
We call expressions of the form 
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the evolutionary vector fields. The expressions ij'/^'^ arc called characteristics of them. 
The value of the evolutionary vector field on a local functional is given by formula 



It is a straightforward calculation to check that this operation is compatible with the 
formal integration by parts, i.e. 

^Div(/) = Div(V/), 

similarly to in the standard formal variational calculus. This relation is, of course, valid 
for integrands. 

It is easy to check that the evolutionary vector field with coefficients 



can be considered as the commutator of the evolutionary vector fields ^ and 77 

VF=K,r?]F = C(r?F)-7?(eF), 

with the Jacobi identity fulfilled for the commutator operation. Therefore the vector fields 
form a Lie algebra. 

^.2. Differentials and functional forms. The differential of a local functional is simply 
the first variation of it 

dfW 



here and below 5(j)^^^ = DjiScpA- It can also be expressed through the Prechet derivative 
or through the higher Eulerian operators 

dF = ^l e^^^ f<^y {5<t>) = Y,j 6^'^DK{Ef{f^'))5<t>A). 

This differential is a special example of functional 1-form. A general functional 1-form 
can be written as 



{J) gAK) 
ak°'Pa ■ 



Of course, the coefficients are not unique since wc can do formal integration by parts. 
Let us call the following expression the canonical form of functional 1-form 



Analogously, we can define functional m-forms as integrals, or equivalence classes 
modulo formal divergences, of vertical forms 



Define the pairing of an evolutionary vector field and 1-form as 

(1) a{0=Uo. = Y.j(^''^'''-AKDKi'I>. 
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The interior product of an evolutionary vector field and functional m-form will be 
^ J " = i Ei-^y^' I ^''^''<^'iU,...A^K^D^ATMl^ A . . . 

Then the value of m-form on the m evolutionary vector fields will be defined by formula 

It can be checked by straightforward calculation that 

Div(a)(a, . ..,U) = Div(a(a, . • 
The differential of m-form given as 

satisfies standard properties 

d' =0 

and 

da(6, • • • , U+i) = Y.^-^y^HM^i, • . . . • . ,^m+l)+ 

i 

+ E(-i)'+'a(fe,^j], a, . . . , ^i, . . . ,4, . . 

i<j 

The Lie derivative of a functional form a along an evolutionary vector field ^ can be 
introduced by the standard formula 

L^a = CJda + d(^ J a). 



3.3. Graded differential operators and their adjoints. We call linear difi'erential opera- 
tors of the form 

N,„ 

'ab 

J=0 Ar=0 

graded differential operators. 

Let us call linear differential operator /* adjoint to / if for arbitrary set of smooth 
functions /a, Qa 



E / fAlABQB = E / SaIXb/b- 



For coefficients of the adjoint operator we can derive the expression 

Krrta^ min{K,J) 

^K-L-M^BA 



K=0 L=0 \ / \ / 

It is easy to check that the relation 

I{x)5{x,y)=t{y)6{:x,y) 

is valid. For example, we have 
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Operators satisfying relation 

i* = -/ 

will be called skew-adjoint. With the help of them it is possible to express 2-fornis (and 
also 2-vectors to be defined below) in the canonical form 



2« 



f>A A IabS4>b- 



It is clear that wc can consider these representations of functional forms as formal decom- 
positions over the basis derived as result of the tensor product of 64) a, with the totally 
antisymmetric multilinear operators 



being coefficients of these decompositions. 

3.4. Multi-vectors and Schouten-Nijenhuis bracket. Let us introduce the dual basis to 
{^(Pa) by formal relation 

and construct by means of the tensor product a basis 

5 5 5 



HbAv) HbAv) "' Hs^iy)' 
Then by using totally antisymmetric multilinear operators we can define functional 
m-vectors (or multi-vectors) 

* = ,„,„i,,.jA-A...AD.„^. 

Here a natural question arises: what is the relation between evolutionary vector fields and 
1-vectors? Evidently, evolutionary vector fields lose their form when integrated by parts 
whereas 1-vcctors conserve it. It is possible to prove the following Proposition [Sol94]. 

Proposition 1. There is a one-to-one correspondence between evolutionary vector 
fields and functional 1-vectors. The coefficients of 1-vector in the canonical form ^^''^ are 
equal to the characteristic of the evolutionary vector field. 

It is not difficult to show that we can define pairing (interior product) of 1-forms and 
1-vectors and this pairing preserves the identification 



When 1-vector is in the canonical form this result coincides with Eq.(l). 

The interior product of 1-vector and m-form or, analogously, of 1-form and m-vector 
is defined as 



^^AU,...,A,nKj<l>T:^^--- 
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Then wc also can define the value of m-form on m 1-vectors (or, analogously, m-vector 
on m 1-forms) 

where in this trace each entry of multilinear operator a acts only to the one corresponding 
^, whereas each derivation of the operator ^ acts on the product of a and all the rest of 

It is possible to extend the differential onto m-vectors 

and analogously onto mixed objects. Evidently, d^tl> = 0. 

With the help of the previous constructions we can define the Schouten-Nijenhuis 
bracket as follows 

for two multi-vectors of orders p and q. The result of this operation is p + q — 1-vector 
and it is analogous to the Schouten-Nijenhuis bracket in tensor analysis [Nij]. Its use in 
formal variational calculus is described in [Dorf]. However this bracket is defined only for 
operators there. We can recommend [01v84] as an interesting source for treatment of the 
Schouten-Nijenhuis bracket for functional multi- vectors. Our construction of this bracket 
guarantees compatibility with the equivalence modulo divergences 

[Div(C),^]5^ = Div[e,ry]^^ = [C, Div(^)] 

Proposition 2. The Schouten-Nijenhuis bracket of functional 1-vectors up to a sign 
coincides with the commutator of corresponding evolutionary vector fields. 



Proof. Let us take two 1-vectors in canonical form 

'Vb 



and compute 

[l^,v]sN = <iUrl-dvJC 

We have 



OiJ) 



and 



Therefore, we obtain 

Proposition 3. (Olver's Lemma [01v86]) The Schouten-Nijenhuis bracket of the 
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two bivectors can be expressed in the form 

where the two differential operators I, K are the coefficients of the bivectors in their 
canonical form. 

Proof. Let us consider the Schouten-Nijenhuis bracket for the two bivectors and 
without loss of generaUty take them in the canonical form 

where = S/ScpA and operators / , K are skew-adjoint. Then we have 

J ocp^ 

and 

-iE / ^(^+^^I)p(^e^ AD..e.j AD.(Cc4T)- 
Now let us make integration by parts in the second term 

dxJV' = -jE / ^^^+^)aA(/ir)'(^<^^c) AD^^s- 

--Jlj 0''^''^''\-^r'{^Q)^-^UAD,^sADj^,_Q{^^^^^^^ 

At last we change the order of multipliers under wedge product in the second term, make 
a replacement M — > M — Q and organize the whole expression in the form 

dX J ^ = - J E / ^^^^""^^^ A {lfl^)'c {K^^,^^n+ 

+(-1)1^1 Q ('^^'^)^P-0-«4"i-^^"^HCc) ADr^is. 

Having in mind the definition of adjoint operator (2) we can represent the final result of 
the calculation as follows, 

[^^^]SN = -^E/ A + AC, 

therefore supporting in this extended formulation the method, proposed in [01v86] for 
testing the Jacobi identity. 
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S.5. Poisson brackets and Hamiltonian vector fields. Let us call a bivector 

formed with the help of the graded skew-adjoint differential operator 

the Poisson bivector if 

The operator I^b is then called the Hamiltonian operator. 

We may call the value of the Poisson bivector on differentials of the two functional 

F,G 

{F,G} = *(dF,dG) = dGjdFj^- 
the Poisson bracket of these functional. 

The explicit form of the Poisson brackets can easily be obtained. It depends on the 
explicit form of the differential of the functional, which can be changed by partial in- 
tegration. Of course, all the possible forms are equivalent. Taking the extreme cases we 
have the expression through Frechet derivatives 

(3) {F,G} = ^le('^Tr{fJi'^g'B), 
or through higher Eulerian operators 

(4) {F,G} = ^J e('^Dp^Q{E^if)ii'^E%)). 

Theorem 1. The Poisson bracket defined above satisfy Definition 1. 

Proof. It follows from the three facts: l)from the previous formulas (3), (4) it is 
clear that {F^ G} is a local functional, 2) antisymmetry of {F, G} as a consequence of 
skcw-adjointness of Iab and 3)equivalence of the Jacobi identity to the Poisson bivector 
property can be proved [Sol94]. 

The result of interior product of the differential of a local functional H and the Poisson 
bivector (up to the sign) will be called the Hamiltonian vector field (or the Hamiltonian 
1-vector) 

idH = -dHJ^ 
corresponding to the Hamiltonian H. 

Evidently, the standard relations take place 

{F, H} = d F{i dH) = (/ d H)F 

Theorem 2. The Hamiltonian vector field corresponding to Poisson bracket of the 
Junctionals F and H coincides up to the sign with commutator of the Hamiltonian vector 
fields corresponding to these Junctionals. 

Proof. Consider the value of commutator of Hamiltonian vector fields IdF and 
I dH on arbitrary functional G 

[/dF,/di/]G = /di^(/di?(G)) -/di/(/dF(G)) = 
= J d F{{G, H})-id H{{G, F}) = {{G, H}, F} - {{G, F}, H] = 
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= -{G, {F,H}} = -Id{F,H}{G), 

where we have used the Jacobi identity and antisymmetry of Poisson bracket. Due to 

arbitrariness of G the proof is completed. 

Example 1. Let us consider a first structure 

{u{x),u{y)} = DJ{x,y) 
of the Korteweg-de Vries equation ([01v86] Example 7.6) 

Uf — '^xxx * 

Construct the adjoint graded operator to 6D according to Eq.(2) 

{OD)* = -6D - D6, 

and the skew-adjoint operator is 

J = i (oD - {eD)* \ =eD + ^De. 



The Poisson bivector has the form 



f 



— AD— 
6u 5u I 



Differential of a local functional H (for simplicity we suppose it is written in canonical 

H = j eh 



form) is equal to 



dH = J eh'{5u) = £ / e^''\-lfE''{h)5u, 

fe=0 



where Frcchct derivative or higher Eulcrian operators can be used. Therefore, the Hamil- 
tonian vector field generated by H is 



or 



or also 



E''{h)D - DE^{h) 



E''{h)D-DE^{h) 



5_ 

d 



The value of this vector field on another functional F coincides with the Poisson bracket 
-dFJdHj^ = {F,H}^^J2j ^k+i (E\f)DE\h) - E\h)DE\f)y 



4. Dynamics on a boundary: an example. Let us now consider a simple Hamil- 
tonian system in order to obtain boundary equations with the help of the new brackets. 
The nonlinear Schrodinger equation 

= --0" + 2ki}\ip\^, 
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can be treated [LR] as generated by the Hamiltonian 

H =^ J{n + n)dx, 



where 



n = r'q' + kr^q^, H = r'q' + kf'^q^ , 



and Poisson brackets are 

{q{x),r{y)} = -2i5{x, y), {q{x), r{y)} = 2iS{x, y). 
To return to the standard form of this equation we should put reahty conditions 

ij) = q = f, ^ = r = q. 
Let us calculate the full variational derivatives 

analogous formulas take place for the bar variables. 

The natural boundary condition arises if we put (5-function contribution on the bound- 
ary to zero by taking 

q' = r' = q' ^f' = 0. 

By considering Poisson brackets for integrals of the total spatial derivatives of canonical 
variables (j)A = {q, r) with the Hamiltonian we expect to obtain dynamical equations on 
the boundary in functional form 

1 1 <j)j^dx = {j <j>A'dx,H} = J D{\Y,lAB^)dx. 



Using Newton-Leibnitz formula we get 



2 



dn 



1 2^ 5<A 



B 



2 



This is different from the standard (bulk) equations 

<pA = lj2^ABE%{n). 

B 

For the given Hamiltonian the formal equations for boundary values are 

qb = -2ikrbqb^, n = 2ikrb^qb, 

if we assume independent behaviour at the ends. It is remarkable that the boundary equa- 
tions are different from the bulk ones despite the boundary condition. These equations 
can be easily integrated and give elementary oscillations at the ends 

Mt) = V^^,(0)e,Tp(-2^fc|^b(0)|2^), 

where the initial value of tp determines both amplitude and frequency of the oscillator. 

In this case the dynamics of boundary values is separated from the bulk dynamics. 
Of course, this situation is not general. 

Let us mention that our boundary condition is compatible with the dynamics, for 
example, 

lq'\l^{Jq",H} = -2tkr'\l=0. 
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5. Discussion. There are not so many publications on problems where divergences 
play a nontrivial role in Hamiltonian formalism of field theory. After the classical paper 
by Regge and Teitelboim [RT] we can recommend the work by Jezierski and Kijowski 
[JK] (see also book[KT]) where the main criterion is also the disappearance of surface 
terms in the first variation of the Hamiltonian. Such functionals are called admissible or 
differentiable, but to be convinced in the full consistency of the formalism it is necessary 
to check the two points: 

1) the space of admissible local functionals should be closed under the Poisson bracket; 

2) the Jacobi identity should be fulfilled for arbitrary admissible functionals; 

and this check is not explicitly demonstrated in the cited works. For infinite domain 
the first requirement was studied by Brown and Henneaux [BH]. In the finite domain 
case the important contribution was made by Lewis, Marsden, Montgomery and Ratiu 
[LMMR] who showed that the standard bracket did not fulfil the Jacobi identity and 
proposed to modify it by adding some surface terms. Unfortunately these authors do 
not consider explicitly the first requirement. We hope that our results could serve as a 
generalization of their ansatz and could be useful in dealing with interesting problems of 
field theory. 
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